UNCLASSIFIED 


OPTICAL  GAIN  PHASE  SHIFT  AND  PROFILE  IN  FREE-ELECTRON 
LASERS(U)  NAVAL  RESEARCH  LAB  WASHINGTON  DC 
B  HAFIZI  ET  AL  20  JUL  87  NRL-HR-5998 


£JU.-U£1 


Uim 


Naval  Research  Laboratory 


Washington,  DC  20375-5000 


NRL  Memorandum  Report  5998 


Optical  Gain,  Phase  Shift  and  Profile 
in  Free-Electron  Lasers 


B.  Hafizi 

Science  Applications  International  Corporation 
McLean,  VA  22102 


P.  Spr angle 

Plasma  Theory  Branch 
Plasma  Physics  Division 

A.  Ting 

Berkeley  Research  Assoc.,  Inc. 
Springfield,  VA  22150 


DTIC 


-iLECTE 
AUG  0  3  1987 


I 


July  20,  1987 


Approved  for  public  release;  distribution  unlimited 


& 


SECURITY  Classification  of  t ms  page 


REPORT  DOCUMENTATION  PAGE 


la  REPORT  SECURITY  CLASSIFICATION 

UNCLASSIFIED 


2a  SECURITY  CLASSIFICATION  AUTHORITY 


2b  DECLASSIFICATION  /  DOWNGRADING  SCHEDULE 


4  PERFORMING  ORGANIZATION  REPORT  NUMBER(S) 

NRL  Memorandum  Report  5998 


3  DISTRIBUTION/  AVAILABILITY  OF  REPORT 


Approved  for  public  release;  distribution 
uni imi ted . 


S  MONITORING  ORGANiZAT'ON  REPORT  NUM8ER(S) 


6a.  NAME  OP  PERFORMING  ORGANIZATION  6b  OFFICE  SvMBOL  7a  NAME  OF  MONITORING  ORGANIZATION 
..  ,  ..  .  ,  ,  (if  applicable) 

Naval  Research  Laboratorv 

Lode  4790 


6c  ADDRESS  «Gfy,  State,  and  ZIP  Code) 

Washington,  DC  20375-5000 


7b  ADDRESS  (City.  State,  and  ZIP  Code) 


8a.  NAME  O!  CUND|NG  SPONSORING 
ORGANIZA7  ON  Armv 

Strategic  Deft’ list-  Communication 


8c  ADDRESS  (C/ty  State,  and  ZIP  Code) 

Huntsville,  Al.  33807-380] 


"  : .  L  tinchjde  Security  Classification) 


8b  OFFICE  SvMBOL  9  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 
(If  applicable ) 


'0  SOURCE  OF  FUNDING  NUMBERS 


PROGRAM  I  PROJECT  TASK 

ELEMENT  NO  I  NO  .v  J  ! !’ i'l)-|  NO 

8322 1C  |7-D-*039  I 


WORK  UNIT 
ACCESSION  NO 


Optical  (.a  in ,  Phase  Shiit  and  I’rolile  in  F'ree-lilect  ton  Lasers 


w  personal  author(S) 

Hafizi,*  B. ,  Spran 


1 3a  tvPE  OF  REPORT 

Inter  im 


le,  P.  and  Ting,**  •' 


1 3b  TIME  COVERED 
FROM _  TO 


1 A  DATE  OF  REPORT  {Year  Month.  Day)  IS  PAGE  COUNT 

I  *»8  ;  .lulv  20  31 


16  SUPPLEMENTARY  NOTATION 

*Sc ience  Appl  i cat  ions  Inti.  C..r».  ,  S!<  I  a-.,  VA  22102 
**Berkelev  Kesearcli  Assoc.,  Inc.,  Spring: icld.  VA  22130 


_ CQSATi  CODES _ I  18  SUBJECT  ’F  RMS  Continue  on  reverse  if  necessary  and  identify  by  block  number) 


SUBGROUP 


-<•  i  * -cl  ron 


p  roi  i  1 1 


R. nli.it  ton  i  mvlopt.'  imirpntuss 

'  >  :  i  luh.iv  ier 


'1  AR$\Ra<  r  Continue  on  reverse  if  necessary  and  identify  by  block  number) 

I  he  uam.  phase  shift.  wavelrmit  .urv.jiurc  .ind  radius  ol  the  radiation  envelope  m  a  ■  ree-elec  iron  laser 
amplifier  are  obtained  m  *hc  small  Men. il  >  ime  I  he  electron  beam  in  assumed  lo  have  a  (iaussian  density 
disiribuhon  m  :h:  c.iry.-s.  NumerK.il  ealeulations  indicate  that  the  radius  and  curvature  of  the 

radiation  beam  cmU  '  10.1’  »  wee  La  isvm;  imt  to  upupie.  spatially-constant  values  alter  a  finite  runsiiion  reuton 
However,  m  the  i^v  mt^niii  tecum  the  w  i-. ..  trouts  are  divergent  \nalytical  exprcismiK  for  the  gam.  phase 
shiit.  >.  u  t  s  a '  ’  1 1  .■  ml  soul  M/e  are  derived  it  :s  shown  analytically  that  small  perturbations  of  the  optical  waist 
ami  wurvatute  a'p»*i:i  the  matched  value  are  spatially  damped  out.  indicating  the  stahihtv  ol  the  matched 
envelope  When  the  eicstton  beam  uivo  n|K  :s  modulated  m  space,  the  optual  spot  M/e  osullales  with  an 
almost  ■  len'i.  .  v\avc.enuih  but  >  -te.r.  eu  o  phase  In  the  uise  ol  small  amplitude  lone  waveienuth  betatron 
im «du la !  1 1 n  the  e!uir>'M  beam  vn,  op.,  ue'ier  ition  ol  optkal  sidebands  m  wave  number  space  is  examined 
and  Lie  ,M:i!  *m  she  if  sv:  m.  •  It  e  *  □  "t  k  s  •  »l  the  pmuarv  wave  is  found  u»  he  neuhuihle  tor  ivpual  evpen* 


D  -  A'.  Mi  TED 


Oc  ABS'RAC^ 

□  /.V-  A, 


NAME  OF  yp*;.-'- 


DO  FORM  1473,V.  va* 


21  ABSiRACT  $ECURI'rv  CL A$S'FiCATiQ\ 

□  otic  USERS  I'NCLASS  1  F 1  CD 


22r>  Tf if puONE  (include  Area  Code)  22 c  <D F F 'C b  SYMBOi 
(2HJI  .  ...  <  <  "do  .;«<) 


S3  A^R  prfi^io”  "■'ay  he  uf.('o  -i"t  i  *  ^"dusted 
A  l  t?d  ?  <*re  :  voImV 


SECoR'lY  CLASS'h'CATiON  of  PAGF 


CONTENTS 


INTRODUCTION  . 

MATHEMATICAL  FORMULATION 

NUMERICAL  RESULTS  . 

ANALYSIS  OF  RESULTS  . 

CONCLUSION  . 

ACKNOWLEDGMENT  . 

REFERENCES  . 

APPENDIX  A  —  Source  Term  . 

APPENDIX  B  —  Sideband  Generation 
APPENDIX  C  —  Stability  Analysis  .... 


OPTICAL  GAIN,  PHASE  SHIFT  AND  PROFILE 
IN  FREE-ELECTRON  LASERS 


I.  Introduction 

A  well-known  feature  of  the  free-electron  laser  (FEL)  is  that  the 

refractive  index  of  the  medium  is  a  complex  function  and  hence  the 

radiation  is  amplified  and  to  some  extent  focused  in  the  vicinity  of  the 
1  2 

electron  beam.  '  It  may  then  be  possible  for  the  electron  and  radiation 
beams  to  interact  over  an  extended  length  along  the  wiggler,  with  the 
diffractive  tendency  being  compensated  by  the  FEL  interaction,  thereby 
enhancing  the  efficiency  of  the  process. 

Considerable  progress  has  been  made  in  studying  this  process  by 
3-8 

several  authors.  The  purpose  of  this  paper  is  to  apply  the  formalism  of 
the  Gaussian-Laguerre  modal  source  dependent  expansion  (SDE)  of  Ref.  8  to 
examine  the  propagation  and  guiding  of  the  optical  wave  in  an  amplifier 
operating  in  the  exponential  gain  regime,  for  a  variety  of  operating 
conditions. 

The  plan  of  this  paper  is  as  follows.  In  Section  II  the  formalism  of 
the  SDE  is  employed  to  obtain  the  evolution  equations  for  the  radius  and 
the  curvature  for  the  lowest  order  mode  of  the  optical  beam,  along  with  the 
relevant  dispersion  relation  for  a  Gaussian  electron  beam  driving  an  FEL 
amplifier  in  the  small  signal  regime.  In  Section  III  numerical  solutions 
of  the  single-mode  equation  for  the  radius  of  the  optical  beam  are 
presented  and  compared  to  the  result  from  a  multi-mode  truncation  of  the 
radiation  field.  In  this  case,  and  for  cases  not  presented  herein,  the 
single-mode  and  multi-mode  results  indicate  that  the  radiation  beam  profile 
entering  the  wiggler  asymptotes  to  a  unique  form  after  an  initial 


transient.  Additionally,  the  numerical  values  of  the  radius  of  the 
radiation  envelope  and  of  the  wavefront  curvature  are  in  fair  agreement, 
irrespective  of  the  degree  of  mode  truncation,  indicating  the  usefulness  of 
the  single-mode  equations.  Limiting  ourselves  to  these  equations,  the 
electron  beam  is  then  allowed  to  oscillate  at  the  betatron  wavelength  and 
the  resulting  radiation  profile  examined.  It  is  found  that  the  optical 
beam  envelope  follows  that  of  the  electrons  with  almost  identical 
wavelength,  but  retarded  in  phase.  Section  IV  discusses  the  results, 
deriving  formulae  for  the  matched  radiation  beam  profile  (i.e.,  radius  and 
curvature)  in  terms  of  the  electron  beam  and  wiggler  parameters.  It  is 
shown  analytically  that  perturbations  of  the  profile  are  spatially  damped 
out,  consistent  with  the  numerical  observations  indicating  a  unique, 
asymptotic  matched  radius  and  curvature.  Appendix  A  presents  the  necessary 
details  required  to  derive  the  source  term,  for  the  wave  equation,  for  a 
planar  wiggler  and  an  electron  beam  with  uniform  density  along  the 
direction  of  propagation.  Appendix  B  considers  the  effect  of  the 
modulation  of  the  electron  beam  on  the  optical  wave.  Specifically,  a 
simple  analysis,  taking  into  account  sideband  generation,  indicates  that 
the  dispersion  characteristics  of  the  primary  wave  are  only  slightly 
modified  for  typical  experimental  parameters.  Appendix  C  presents  the 
details  of  the  stability  calculation. 


II.  Mathematical  Formulation 

The  purpose  of  the  present  section  is  to  present  the  salient  features 

g 

of  the  source  dependent  expansion  method  so  as  to  fix  the  notation  and  for 
reference  in  the  subsequent  sections. 

For  a  planar  wiggler,  it  is  appropriate  to  assume  a  linearly  polarized 
radiation  vector  potential 


with  angular  frequency  <o  and  complex  amplitude  A.  In  the  slowly  varying 
envelope  approximation,  the  wave  equation  reduces  to 


13  3  1  3^  2iw  ^  > 

r  3r  r  3r  +  ~2  772  +  “c  3z  a  =  s(r»e»z>» 
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where  a  =  |e|A/mQc  ,  and  the  source  function  is  given  by 
S(r,9,z)  =  -  — -L®  I-  |jx(r,6,z)exp  -i  —  -  w  tj  j 


(1) 


(2) 


Here  e  is  the  charge  on  an  electron  of  (rest)  mass  mQ,  Jx(r,0,z)  is  the 
current  density  and  {  }glow  indicates  that  only  the  spatially  and 
temporally  slow  part  of  the  quantity  in  braces  is  to  be  retained. 

The  basic  premise  of  the  work  presented  herein  is  that  the  radiation 
field  is  azimuthally  symmetric  and  the  vector  potential  is  expressible  as 


00 

a  ( r ,  0,  z)  =  £  am(z)Dma’z)’ 
m=o 


(3) 


with  Dm  =  LJS)  exp  {-[l-ia(z) ] £/2} ,  where  £;  =  2r2/r2(z),  r  (z)  is  related 
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to  the  radiation  spot  size,  a(z)  is  proportional  to  the  curvature  of  the 
wavefront,  and  Lm(£)  is  the  Laguerre  polynomial  of  order  m. 

Now,  if  the  transverse  prclile  of  the  radiation  beam  is  close  to  a 
Gaussian,  the  lowest  order  mode  is  expected  to  dominate'*’'^’®,  and, 

following  Ref.  8,  it  is  simple  to  show  that  the  associated  vector  potential 
evolves  according  to 

(fe  *  \k  •  -1  F0’ 

and  the  spot  size  and  wavefront  curvature  evolve  via 


(4) 


v.fviuivmMr 
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dz 


a  -  2(l+a2) — — 


co  r 


(5a) 


(5b) 


where 


*  1  d  .  r,,  2.  c  ad  Id  1 

Ao  “  r  dz  rs  +  1 L  +a  ^  2  r  dz  rs  +  2  dz  aJ ’ 

s  cor  s 


the  F's  are  given  by  the  following  overlap  integral! 


Fm(z)  "  h  J  S(5,z)Dm*(t,z),  (6) 

o 

and  the  label  R  (I)  indicates  the  real  (imaginary)  part. 

Noting  that  Lo(£)  =  1,  the  normalized  vector  potential  is  seen  to  be 
given  by  (Eq.  (3)] 


a( r , 0, z)  =  a  (z)  exp(-[ l-ia(z) — }, 

l  rf(z)J 


where,  in  the  exponential  gain,  small-signal  regime, 


( ( z )  =  a(o)  exp  |i  JdZj ^Ak(z1)-ir(z1)j j. 


(7) 


(8) 


Here  a(o)  is  the  input  signal  at  z  =  o,  and  the  two  components  of  the 
refractive  index  are  given  by 


,  i  c  Ak ,  .c(t-  29  l-i 

n  =  [1+ -  -  l  -  i  -  r  t  — s— 

z  1  w  J  (0  V  9z  2  J 


(9a) 


2cr  ,  . . 

nr  =  — 2  +  1>* 

cor 
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(9b) 


Assuming  the  electron  beam  profile  to  be  given  by 


n.  (z)  =  n. 
b'  '  bo  r^z) 


rb(z> 


where  r^(z)  is  the  electron  beam  radius  at  z  and  n^  is  the  beam  density  at 
rfe(z)  =  rbQ,  the  source  term  in  Eq.  (1)  may  be  readily  evaluated  (Appendix 
A),  to  obtain 


S(r,z)  = 


p2  “bo  rbo 

B  2y3c2  lrb(2) 


2  'i  w  k  a  a 
r  w  w 


r2  J  c(Ak-ir)2' 


where  the  vector  potential  of  the  planar  wiggler  of  periodicity  2n/kv  is 
given  by 

A  w  =  Awcos(kwz)e  x,  (12) 

aw  -  le  |Aw/™0c2,  (13) 

y  is  the  relativistic  mass  factor,  f0  is  the  usual  difference  of  Bessel 
functions,  =  JQ(0  -  C  =  (l/4)a2/[l  +  (l/2)a2],  and 

(  2  )1/2 
“bo  =  Nel  nbo/moJ 

is  the  plasma  frequency  of  the  electron  beam  with  density  n,  . 

bo 

Substituting  Eqs.  (8)  and  (11)  into  Eq.  (6)  and  making  use  of  Eqs.  (4) 
and  (5),  it  is  simple  to  show  that  the  equations  reduce  to 

<ua, 


4^1' 

d  (kwz) 


ck  )  (F  .  )  .  ,2 

— “  "  2  j - -  fk  r  1  , 

.  “  J  rwa0ji  {v  s) 


(14b) 


\  -r .  s ,  -r,  »,  •  «  v  '  -  "  «  *  **  »  ■ *  ././  *  *  *  4  %’  tAa*.'  *, 


V  f+j '  •*  v 


•  *■  .V  .v.j 


0, 


(14c) 


where 


k  a 
v  o 


(14d) 


The  spatial  evolution  of  the  system  is  governed  by  the  differential 
system  (14a)  and  (14b)  along  with  the  dispersion  relation  (14c),  the 
solution  of  which  yields  a(z),  rg(z),  Ak(z)  and  I"(z). 


III.  Numerical  Results 

Having  obtained  the  single  mode  system  of  Eqs.  (14),  it  is  of  interest 
to  determine  the  extent  to  which  it  approximates  the  general  solution  in 
(3).  Once  it  is  established  that  Eqs.  (14)  provide  an  adequate 
representation  of  the  general  solution,  it  is  then  possible  to  study  a 
variety  of  problems  of  interest  by  solving  a  simple  set  of  equations. 

Briefly,  the  numerical  procedure  for  solving  an  initial-value  problem 
is  the  following.  Substituting  Eq.  (14d)  into  Eq.  (14c)  yields  a  cubic 
(algebraic)  equation  for  &k  -  iT  which  may  be  solved,  at  each  z,  in  terms 
of  rs(z),  a(z)  and  r^(z),  thus  enabling  Eqs.  (14a)  and  (14b)  to  be  stepped 
forward  in  z.  Since  in  the  absence  of  source  terms  an  input  radiation 
signal  diffracts  away  on  the  scale  length  defined  by  the  Rayleigh  range  z„, 


z 


R 


2.  . 

«  rs(z> 

— 2c 


z=0’ 


(15) 


it  is  informative  to  present  the  numerical  results  with  the  distance  along 
the  wiggler  measured  in  units  of  the  Rayleigh  range.  In  all  the  numerical 
results  to  be  presented,  the  radiation  field  is  assumed  to  be  in  the  form 
of  plane  waves  at  the  entrance  to  the  wiggler,  i.e.,  a(z=0)  =  0. 

Case  I 

To  begin  with,  Fig.  1  shows  the  results  for  the  following  parameters: 

beam  current,  I,  =  270  A,  r,  =  0.01  cm,  y  =  2000,  2ti/k  =  10  cm,  a  «  6.15 
b  ’bo  w  w 

1/2 

and  rs(z=0)  =  0.02  cm.  Noting  the  factor  of  2  difference  between  the 
definition  of  ay  in  Eq.  (13)  and  that  in  Ref.  4,  it  is  clear  from  Fig.  1(a) 
that  after  a  transient  oscillation  over  a  distance  of  about  20  Rayleigh 
ranges,  the  radiation  spot  size  approaches  a  value  quite  close  to  that 
obtained  with  the  two-dimensional  FEL  code  FRED  at  the  Lawrence  Livermore 

4 

National  Laboratory  (LLNL).  We  also  find  that  for  all  the  numerical  cases 
examined,  a  unique,  asymptotic  spot  size  is  obtained  irrespective  of  the 
initial  optical  waist.  Figure  1(b)  shows  the  spatial  evolution  of  a, 
indicating  that  it,  too,  approaches  a  constant  value  after  an  initial 
transient  behavior. 

The  solid  curve  in  Fig.  2  shows  the  evolution  of  1/e  width  of  the 
radiation  amplitude  with  a  five  mode  (m=0, 1 , 2 , 3 , 4)  source  dependent 
expansion  calculation  using  the  same  set  of  FEL  parameters.  The  radiation 
field  is  represented  by  Eq.  (3)  and  the  source  term  is  given  by  Eq.  (11). 
With  the  assumption  that  the  fundamental  mode  dominates,  only  the  Ak  and  T 
of  aQ  (r,z)  are  involved  in  the  source  function  and  they  are  obtained  from 
Eqs.  (14c)  and  (14d).  It  is  found  that  the  fundamental  mode  remains 
dominant  over  many  Rayleigh  lengths.  For  comparison  the  dashed  curve  in 
Fig.  2  shows  the  fundamental  mode  spot  size  of  Fig.  1(a),  and  the 
asymptotic  results  are  seen  to  differ  by  about  10%.  This  suggests  that  the 


single-mode  system  of  Eqs.  (14)  may  be  regarded  as  a  reasonably  accurate 
simplification  of  Eq.  (3).  Henceforth  the  results  presented  pertain  to 
Eqs.  (14). 


Case  II 

Figure  3  presents  the  results  for  a  case  where  the  electron  beam  is 
not  matched;  i.e.,  the  envelope  of  the  electron  beam  is  modulated: 


rb(z)  =  rbo  +  Srbsin(k0z)» 


(16) 


where  Sr,  is  the  amplitude  of  the  modulation  and  for  simplicity  ka  is 
D  p 

9  _ 

chosen  to  be  equal  to  the  betatron  wave  number  kvav/(-l2Y0z) ,  neglecting 
self-fields.^  is  the  beam  speed  along  the  wiggler  axis  normalized 
to  c.  The  parameters,  typical  of  the  Advanced  Test  Accelerator  experiment 


at  LLNL,  are  I,  =  2  kA,  r,  =  0.3  cm,  y  =  100,  2n/k  =  8  cm,  a  =  1.72, 
b  bo  w  w 

rs(z=0)  =  0.35  cm.  (The  reader  is  referred  to  Refs.  9  and  11  for  details.) 
In  Fig.  3,  where  Sr^/r  bo  =  O’1’  is  observed  that  the  optical  spot  size 
follows  the  modulations  in  the  electron  envelope  apparently  identically. 
Specifically,  a  number  of  cases  were  examined  with  Sr^/r^  UP  to  0.4.  In 
all  cases  the  electron  and  optical  beams  oscillate  with  almost  identical 
wavelength,  although  the  radiation  beam  appears  to  lag  behind  in  phase. 


Defining  the  modulation  depth  A  =  [(r)  -  (r)  .  l/[(r)  +  (r)  .  1, 

it  is  found  from  Fig.  3(a)  that  Ag  =  0.087  whereas,  from  Eq.  (16), 

A^  =  6r^/r^o  =  0.1.  Although  .he  modulation  depth  of  the  electron  beam 

differs  from  that  of  the  radiation  beam,  it  is  found  that  A  increases 

s 

with  6r,  . 
b 

More  generally,  allowing  for  the  defocusing  effect  of  self-fields, 


there  is  always  the  possibility  of  a  small  amplitude  ripple  on  the  electron 
beam  envelope  and  hence  on  the  radiation  beam  envelope.  In  Appendix  B, 


generation  of  sidebands  is  considered  in  a  simplified  model  and  found  to 
have,  for  typical  cases,  an  insignificant  effect  on  the  linear  dispersion 
characteristics  of  the  primary  optical  wave,  as  implicitly  assumed  by 
employing  the  source  term  in  Eq.  (11)  in  the  present  case. 


IV.  Analysis  of  Results 

One  interesting  feature  of  the  numerical  results  is  that  in  all  cases 
the  radiation  spot  size  has  a  unique,  asymptotic  limit  irrespective  of  the 
initial  value.  The  asymptotic  value  of  rg  and  of  a  is  determined  by  the 
fixed  points  of  Eqs.  (14a)  and  (14b);  i.e.,  at  the  fixed  point 


(17a) 


(17b) 


Combining  Eqs.  (17a)  and  (17b)  one  obtains 

w^o' 


2  ckv 

(1  -  ia;  —  + 


<Vs>2(jrf)  ■ 


which,  upon  making  use  of  Eq.  (14d),  yields 


k2rhn1/2 

*  -  -nn—  r2?’  r ' 4k/a' 

1  +  a 


where 


2 

and  f  =  (r^/rs)  is  the  filling  factor.  Substituting  the  expressions  for 
Ak  and  f  into  the  dispersion  relation  (14c),  one  obtains 
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a  =  [ f / ( 3  f  +  2)] 


^  =  23/4k  yfR1/2 
w  B 


,/v)1M  (U4/2)3/4  f1/4(L2f)3/2 


1/2  "1/4  ’ 

(l+3f/2r/4 


where  v  =  (“bo^c/^c/  is  Budker's  parameter.  These  expressions  may  be 


used  to  obtain  the  asymptotic  spot  size  for  a  given  filling  factor,  and 


then  one  obtains  the  corresponding  electron  beam  radius  via  r^  =  rsfA  . 


To  avoid  complications  arising  at  the  outer  edges  of  the  optical  beam, 


where  the  field  amplitude  is  small,  typically  a  filling  factor  f  <  1/2  is 


appropriate.  It  is  also  possible  to  rearrange  the  expression  for  rs  to 


obtain 


f!  .  £2  .  Q--  |,)f  -  ,  .  0, 


where 


f  2  /Yrbkw)4l1/3  1 

K  Ir/vJl  2  J  J  2  ' 

v 


The  cubic  equation  for  f  may  be  solved  to  obtain  an  explicit  expression  for 


r  .  Noting  that  the  sum  and  the  product  of  the  three  roots  of  the  cubic 


equal  -1  and  q,  respectively,  it  follows  that  there  is  a  unique,  real  value 


for  the  asymptotic  spot  size  rg . 


To  examine  stability,  it  is  convenient  to  define 


.  r 

1  iT’ 

w 


and  substitute  Eq .  (14d)  into  Eq.  (14c)  to  obtain  the  local  dispersion 


rela  t ion : 


10 


Y3  +  2 


(ck 

1-ia  „2 

ck 

(r.  ^ 

w 

2  _ —  p 

b 

(A) 

2  1 

(k  r  )Z 
'  w  s 

2  <*>  n 

r 

s) 

I+CVV 


-  •  "I'.J  KvT' 

which  may  be  solved  iteratively.  It  turns  out  that  for  the  parameters  of 
Case  I,  at  the  lowest  order,  the  right-hand  side  balances  the  quadratic 
term  on  the  left.  The  relevant  root,  with  Ak,  T  >  0,  may  be  substituted 
into  Eq.  (14d)  to  obtain,  for  a  >  o, 


-ia _  1  f  n 

2  2  2  2 


(k  r  )  +(k  r,  ) 
'  w  s  v  w  b' 


■  ,1/2  a-i  l+(l+a2}/2 

2(l+a2)  l+(l+a2)1/2  1/2 


(k  r  ) 

v  w  s  _ 

(k  r  )2+2(k  r,  )2 
v  w  s  v  w  b 


(k  r  ) (k  r.  ) 
v  w  s'  w  b 


wrb)  f(k  r  )2+(k  r.  )21 
Lvws  vwb'J 


Perturbing  Eqs.  (14a)  and  (14b)  about  the  fixed  point  and  making  use  of 
Eq.  (19),  it  is  simple  to  show  that  the  perturbation  is  spatially  damped, 
thus  indicating  the  stability  of  the  fixed  point.  The  algebraic  details 
are  relegated  to  Appendix  C. 

Another  aspect  of  the  results  which  is  of  interest  pertains  to  the 

nature  of  the  phase  fronts  and  the  flux  of  optical  power  in  the  asymptotic 

region.  From  Eqs.  (7)  and  (8)  it  is  simple  to  check  that,  in  differential 

2 

form,  the  surfaces  of  constant  phase  are  given  by  (co/c  +  Ak)$z  +  (2ra/rs)5r 

=  0,  and  hence,  noting  that  Ak,  a  >  0,  the  wavefronts  are  divergent  in  the 

direction  of  propagation.  Consistent  with  this,  there  is  a  nonvanishing 

transverse  component  of  the  Poynting  flux.  Specifically,  for  <  1  the 

ratio  of  flux  of  optical  energy  in  the  transverse  direction  to  that  along 

2 

the  z  axis  is  ~  ar/kr  <<  1. 

s 


V. 


Conclusion 


Based  on  the  results  presented  herein,  the  simplicity  and  accuracy  of 
the  single-mode  Gaussian-Laguerre  approximation  to  the  solution  of 
Maxwell's  equations  has  been  demonstrated.  It  is  shown  that,  in  the 
exponential  gain  regime  of  operation  of  an  FEL  amplifier,  there  is  a 
unique,  asymptotic  spot  size  for  the  radiation  beam  irrespective  of  that  at 
the  entrance  of  the  wiggler.  There  is,  however,  a  transverse  flux  of 
optical  power.  It  is  shown  analytically  that  the  asymptotic  profile 
(i.e.,  the  radius  and  the  curvature  at  large  z)  is  stable  to  small 
amplitude  perturbations.  With  a  spatially  modulated  electron  beam 
envelope,  that  of  the  optical  beam  is  found  to  oscillate  on  the  same 
spatial  scale. 
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Appendix  A:  Source  Term 


In  this  appendix,  the  details  of  the  evaluation  of  the  source  term  S 
in  Eq.  (11)  are  presented. 

The  FEL  source  current,  Jx(r,0,z),  in  a  linear  viggler  is  given  by 


Jx(r,0,z)  =  -  |e  j&n^r.ejZ)^ 


i  i2t  -ik  z 

-  Ie I  ^nbe  w 

2vni  c 
o 


A  +  c .  c . 
v 


where  6n^  is  the  perturbed  beam  density  and  the  relation  v^  =  vy  =  |e|Ay  * 
cos(kvz)/>moc  has  been  used.  Equation  (2)  can  then  be  written  as 


S( r , 0, z) 


'4n|e|2Snbaw  -i  [<k  +  ky)z 

- 2 -  e 

'  ^0° 


(Al) 


where  k  =  w/c. 

The  perturbed  beam  density  can  be  evaluated  from  the  continuity 
equation, 


d6n,  36v 

_ b  _  _ z 

dt  “  nb  3z  ’ 


(A2) 


and  the  equation  of  motion  in  the  z-direction, 


(A3) 


where  electron  self-field  effects  are  neglected.  Taking  the  convective 
time  derivative  of  Eq.  (A2),  and  incorporating  the  linearized  version  of 
Eq.  (A3),  one  can  arrive  at  the  following  equation  for  the  perturbed  beam 


Ste 
I*$i 


•IjSs 

/•ill 


yin  3z  (3z  2  3t)  pond 

o  c 


where 

_  lelAvA  i[(k  +  k  )z  -  «tl 
*  .  =  - ~  el  w7  J  +  c.c. 

P°nd  Aym  c2 
o 

With  the  assumption  that  A(r,9,z)  is  a  slowly  varying  function  of  z, 
i.e.,  |31nA/3z|  <<  k  <<  k,  Eq.  (A4)  becomes 

i  y 


d2  &n,  |e  |2n,  A  A  i 

_ b  _  _ _ ! _ b_w_  k  k 

,2  "  .  2  2  2  w  e 

dt  2y  m  c 

o 


1  [(k  +  kw)Z_Wt] 


where  the  resonance  condition,  «  =  v  (k  +  k  )  is  used. 

’  z  v 

For  a  near  Gaussian  radiation  field  in  the  exponential  gain  regime, 


A(r,0,z)=A  (r,9,z)=A  (o)expji  [ Ak(z. )-ir(z.. ) ]dz. - [ l-ia(z) ]  - 

°  °  U  1  t(z) 


and  assuming  Ak,  T,  a  and  rg  are  slowly  varying  functions  of  z,  Eq.  (A5) 
can  be  integrated  immediately  to  give 

|e|2nbAAkk  i[(k  +  k  )z  -«  tl 

&n.  =  — * — - w  e  l  w  j  +  c>c>  6 

b  n  2  24,..  ...z 

2y  m  c  ( Ak-iT) 
o 

When  Eq.  (A6)  is  substituted  into  Eq .  (Al),  taking  into  account  the 
usual  difference  of  Bessel  functions  for  a  planar  wiggler,  and  Eq.  (10)  for 
the  beam  profile,  the  source  function  in  Eq.  (1)  is  then  given  by  Eq .  (11). 


wv 

mnlRR 
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Appendix  B:  Sideband  Generation 


In  this  appendix  generation  of  sidebands  to  the  primary  optical  wave, 

due  to  the  spatial  modulation  of  the  electron  beam,  is  analyzed.  It  is  to 

be  emphasized  that  the  following  analysis  is  intended  merely  to  show  that 

the  dispersion  characteristics  of  the  primary  optical  wave  are  only 

2 

slightly  modified  [~  (SNq/No)  ]  for  typical  experimental  parameters,  as 

implicitly  assumed  in  applying  the  results  of  Appendix  A  to  the  case  of  a 

modulated  electron  beam  in  Section  III. 

The  development  of  the  linear  theory  herein  generalizes  that  of 
12 

Sprangle  et  al.,  to  which  reference  should  be  made  for  further  details. 

The  form  of  the  vector  potential  of  a  planar  wiggler  employed  in  this 
appendix  is  slightly  different  to  that  given  by  Eq.  (12): 


A  =  A  exp(ik  z)  -  c.c.  e  , 
-w  w  L  w  J-x 


where  A  is  purely  imaginary,  and  that  of  the  linearly  polarized  radiation 

V 

field  is  taken  to  be  of  the  form 


A  =  |a+  expj^ik+z-iw  tj  +  A  exp  ^ik_z-iwt j 
+  Aq  exp^ik  z-iw  tj  +  c.c.j  ex, 


where  it  is  assumed  that  the  electron  density,  modulated  at  the  betatron 

wavelength  2n/ka,  has  the  simple  form 
P 


SN 

n  =  N  +  -=■ 

OOl 


-  exp(ikgz)  +  c.c.j, 


with  k0  <<  k  <<  k,  and  k  =  k  +  kOI  k  =  k  -  k„. 
gw  +  g  -  g 
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r>'*, 


■  *■  m  *#•  *  *  *  •  w  m  * 


Following  Ref.  12,  the  wave  equation  is  found  to  be 


(2_  _  L_  iL  _  !^_J|  A  .  inkl 

3z2  c2  3t2  y  c2  y  m  c- 

o  'o  o 


2  &nV 


where  yq  is  the  relativistic  factor  in  the  absence  of  the  radiation  field, 
=  (4nno  |e  |  /mQ)  ,  and  6n  is  the  density  perturbation  caused  by  the 
radiation.  Note  that  the  velocity  vzq  along  the  wiggler  axis  is  not 
affected  by  the  betatron  oscillation  and  hence  yq*  to  lowest  order  in 
|eAw/Yomoc  |  ,  is  not  a  function  of  z.  Defining  the  ponderomotive 


potential 


pond 


del 


Y  m  c 
o  o 


2n  •  n 

-w  - 


A, 


the  momentum,  continuity  and  Poisson's  equations  may  be  combined  to  obtain 


,3n 


d  5n  _  o  d 

dt2  no  Uz  Jdt 


4Ttnn|e|- 


m  y  Y_ 
o  o  z 


Sn  + 


m  y  Y 
o  o  z 


f!d>)  8* 

Uz  J  3z 


_  d lL  <L  n  Vzo  3 

m  y  8z  o  (3z  +  2  3tJ  pond* 

o  o  c 


(Bl) 


2  2  -1/2 

where  yz  =  (1-v  /c  )  ,  ♦  is  the  scalar  potential,  and  terms  such  as 

2  2  2 
3  nQ/3z  ,  which  are  on  the  order  of  k^,  have  been  neglected. 

Writing  .  k  .  ks,  k_  .  k  -  ke, 


&n 


=  j&n  +  exp  i  ^k+  +  k^jz-ioo  tj  +  Sn  exp  i  ^k  +  kvjz-iw  tj 


Sn 


oexp  [i  (k  .  ky)r  -icotj  +  c.c.j, 


noting  that,  on  the  left-hand  side  of  Eq .  (Bl),  the  ratio  of  the  fourth  to 

the  third  term  is  on  the  order  of  k_/k  <<  1,  one  finds  that 

P 
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where  e  =  (6Nq/2No),  and  m^,  m^  =  0(e  ).  It  is  then  simple  to  show  that 
2 

correct  to  0(e  ),  the  dispersion  relation  is  given  by 


Yo 


is  the  usual  matrix  element  for  the  primary  wave, 
=  m22^k-^ ’ 

m12  =  mi2(k+’  k’  kp) 


a 


^bo 


k„v_ 


+  y  c 
o 


2  re  zo 


w 


-  fk  +  k  )v 
l.  v)  zo 


^bo 


YoV 


2  1/2 

and  w,  =  (4ne  N  /m  )  .  Note  that  with  the  definition  chosen  for  A 

bo  o  o 

2  2  2 

this  appendix,  a  =  (eA  /me)  <0. 

v  v  w  o 

To  proceed  along  the  lines  of  Ref.  12,  it  is  convenient  to  write 
m22  =  M22  +  C22 ’ 


where 


■  {[-  -  (*  *  Kho]  -  (*2  -  4  *  % 

v  r0rz)  c  yqc 


‘22 


and 


'22 


o  2 

%o  ,  ,  2 

- k  k  a 

3  w  w 


is  the  "coupling"  term.  The  dispersion  relation  then  becomes 


,  &N  .  ^  -  a  . 

.SN  (a  a 

U  (  2FT  )  +  i^) 

M22  “ 

r  M  (m33  *mll  J 

C22 

•'SV.'C.V 


Mjj  yields  the  dispersion  relation  for  uncoupled  electromagnetic  and  space 
charge  waves.  The  right-hand  side  of  Eq.  (B2)  introduces  the  FEL 
interaction  and  coupling  to  sidebands,  and  its  effect  is  included 
iteratively.  At  the  lowest  order,  M22  =  0  for  some  (w,k).  Substituting  in 
the  right-hand  side,  the  second  set  of  terms  vanishes;  the  term 
proportional  to  C22  survives. 

Substantial  modification  of  this  dispersion  relation  is  expected  if 


m33  mll 


i .e. ,  if 


o  0  ,  -3/2 

“  "  2  kSc“bo 


(2k  vd 


1/2  1/2  3/4 

Y„ 

z  o 


For  typical  experimental  parameters,  the  right-hand  side  of  this  equation 
exceeds  unity,  whereas  6Nq/No  <<  1,  implying  the  insignificance  of  the 
effect  of  modulation  on  the  dispersion  relation. 


Appendix  C:  Stability  Analysis 


The  purpose  of  this  appendix  is  to  establish  the  stability  of  the 
fixed  point  (rg,a)  of  Eqs.  (14). 

Perturbing  Eqs.  (14a)  and  (14b)  about  the  fixed  point  and  making  use 
of  Eq.  (19),  it  is  seen  that  the  perturbation  evolves  according  to: 

d  (6a) 

d(k  2)  UxJ 

W 

-  2  (“11  a  1 2 1  f!«). 


a2i  a22; 


2  2 
where  x  =  (k  r  )  ,  y  =  (k  r,  )  , 
'  w  s  v  w  b 


I  3a 


(XR  '  «*l)- 


-(1+a2) (ckw/«jy(2x+y)  g 


2  2 

x  (x+y) 


I  3  Y 
a21  =  ~a~  ~  X  3a  X!  ’ 


(V  «*i)' 


a  ck  /(*)  -  ,  . 

7^-  kM- 


(x+y) 


(  2i 
Ua  . 

/ 


r  (  2)X/Z1 

a-i  1+  1  +  a  ,  .1/2 

;  2,l,l  ^  fe)  •<«> 

H1**)  J 


Assuming  that  6a,  6x  ~  exp(Xk  z),  one  finds  that 

V 


a (ck^/wj (x+2y) 


(x+y)4 


+  S2  +  + 


■a  ^ckv/wjx 


(x+y) 


2  +  S2  S1 


finds  that  Re  X  <  0,  thus  indicating  the  stability  of  the  fixed  point  to 


small  amplitude  perturbations. 
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